Evolution of the Kondo resonance feature and its relationship to spin-orbit couphng 
across the quantum critical point in Ce2Rhi ^.Co^jSia 
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We investigate the evolution of the electronic structure of Ce2Rhi_3;Coa;Si3 as a function of x 
employing high resolution photoemission spectroscopy. Co substitution at the Rh sites in antifer- 
romagnetic Ce2RhSi3 leads to a transition from an antiferromagnetic system to a Kondo system, 
Ce2CoSi3 via the Quantum Critical Point (QCP). High resolution photoemission spectra reveal dis- 
tinct signature of the Kondo resonance feature (KRF) and its spin orbit split component (SOC) in 
the whole composition range indicating finite Kondo temperature scale at the quantum critical point. 
We observe that the intensity ratio of the Kondo resonance feature and its spin orbit split compo- 
nent, KRF/SOC gradually increases with the decrease in temperature in the strong hybridization 
limit. The scenario gets reversed if the Kondo temperature becomes lower than the magnetic order- 
ing temperature. While finite Kondo temperature within the magnetically ordered phase indicates 
applicability of the spin density wave picture at the approach to QCP, the dominant temperature 
dependence of the spin-orbit coupled feature suggests importance of spin-orbit interactions in this 
regime. 

PACS numbers: 71.27.-l-a, 71.28.-l-d, 75.30.Fv 



I. INTRODUCTION 

Quantum criticality^ in heavy fermion systems is ob- 
tained by tuning the antiferromagnetic phase transition 
temperature to zero by external non-thermal control pa- 
rameter like pressure, chemical composition, magnetic 
field etc. The study of quantum criticality has attracted 
a great deal of attention recently due to the finding of 
non Fermi liquid excitations in the properties of metals 
near the magnetic instability in correlated systems and in 
particular in unconventional superconductors. The quan- 
tum critical point (QCP) separates two entirely differ- 
ent ground states - (i) long range magnetically ordered 
state and (ii) paramagnetic/nonmagnetic Fermi liquid 
state3 While the physics in the ground states on either 
sides, away from the QCP are reasonably well under- 
stood, the physics in the vicinity of QCP is a long stand- 
ing puzzle.—^— The question to be answered is, how the 
quantum critical fluctuations affects the electronic exci- 
tations and how it evolves as one goes from one side of 
QCP to the other side. 

In order to address the issue, we studied an archetyp- 
ical systems Ce2Rhi_2;Co2;Si3 - the compounds in the 
whole composition range crystallize in a AIB2 derived 
hexagonal structure with space group P6/mmm as 
shown in Fig. li^"— The valence band consists of Si 3p, 
Co 3d and Rh Ad electronic statesi^iiS Since the bind- 
ing energy of the 4(1 electronic states is higher than the 
3d electronic states and 3c? electrons are relatively more 
correlated due to smaller radial extension, ^^'^^ the sub- 
stitution of Co at Rh sites influences the effective d-p 
hybridization as well as the itineracy of the conduction 
electrons that leads to a plethora of interesting proper- 
ties. For example, one of the end members, Ce2RhSi3 ex- 
hibit antiferromagnetic ordering at about 7 K.^ Various 




FIG. 1: Crystal structure of Ce2Rhi_a;Rh2,Si3 forming in 
AIB2 type hexagonal structure. 



studies involving lattice compression/expansion via ex- 
ternal pressure or chemical substitutions suggested that 
this compound lies close to the peak of Doniachs mag- 
netic phase diagramJ^ii^ The sample with x = 0.6 com- 
position exhibits quantum critical behavior such as non- 
Fermi liquid behavior in the electrical transport measure- 
ments. The compositions with 0.0 < x < 0.6 exhibit 
antiferromagnetic ordering. It was shown that applica- 
tion of high pressure on the compounds in this compo- 
sition range leads to a change in properties similar to 
that observed via change in composition.^^ The magnetic 
susceptibility measurements on the other end member, 
Ce2CoSi3, does not show magnetic ordering down to 0.5 
K akin to a Kondo lattice compound with an estimated 
Kondo temperature, Tk of about 44 K. The temperature 
dependence of electrical resistivity of Co-end members 
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appear to be similar to a typical mixed valent Kondo 
lattice systemj^siS 

In this article, we present our results of investigation 
on the evolution of electronic structure as a function of 
Rh/Co concentration employing high resolution photoe- 
mission spectroscopy. The high resolution spectra exhibit 
distinct signatures of Kondo resonance feature and its 
spin orbit coupled component in the whole composition 
range. The relative intensity of the Kondo resonance fea- 
ture with respect to the spin orbit split feature exhibits 
distinctly different temperature evolution in the different 
parts of the phase diagram. 



II. EXPERIMENTAL DETAILS 

High quality samples of Ce2Rhi_a;Coa;Si3 for various 
values of x were prepared in the polycrystalline form 
by arc melting - growth of the material from congru- 
ently molten phase ensured large grain size with ho- 
mogenous stoichiometry. The samples are characterized 
by x-ray diffraction, transport and magnetic measure- 
ments as described elsewhere^S' The photoemission mea- 
surements were carried out using high resolution Gamma- 
data Scienta electron analyzer, SES2002 and monochro- 
matic photon sources at a base pressure better than 3 
X 10~^^ torr. The energy resolutions were set to 300 
meV at Al Ka (1486.6 eV) photon energy, and 5 meV 
at He II (40.8 eV) and He i (21.2 eV) photon energies. 
The samples were very hard and could not be cleaved 
or fractured. Hence, the sample surface was cleaned by 
in situ scraping at each temperature using a fine grained 
diamond file and the surface cleanliness was ensured by 
the absence of O Is and C Is signals. The reproducibil- 
ity of the spectra was ensured after each trial of surface 
cleaning. The Fermi level was determined at each tem- 
perature by measuring the Fermi cutoff in the spectral 
function of high purity silver mounted in electrical con- 
tact with the sample. The temperature variation down 
to 5 K was achieved using a open cycle helium cryostat, 
LT-3M from Advanced Research Systems, USA. 
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FIG. 2: (a) XP valence band spectra of Ce2Rhi_a:CoiSi3 for 
various values of a; at 5 K. The solid line shows a typical in- 
tegral background function, (b) Background subtracted high 
resolution valence band spectra exhibiting evolution of the 
spectra with a; at 5 K. (c) The same background subtracted 
spectra shown after normalization at 1.5 eV binding energy to 
demonstrate how the spectral intensity changes with compo- 
sition. The gray shade shows the increase in Co 3d intensity 
and blue shade shows subsequent decrease in Rh 4d intensi- 
ties with the increase in x. LDA results of the partial density 
of states of (d) Ce2RhSi3 and (e) Ce2CoSi3. The t/-scale is 
expanded to visualize the Ce 4/ (yellow shaded region). Si p 
and Co/Rh d partial density of states (PDOS) in the occupied 
part with clarity. The inset in (e) shows the Ce 4/ PDOS for 
Ce2CoSi3 (dashed line) and Ce2RhSi3 (solid line). 



III. RESULTS AND DISCUSSIONS 

The photoemission measurements were carried out on 
four compositions, namely, x = 0.0 (antiferromagnetic) , 
0.3 (30% Co substituted and antiferromagnetic), 0.6 (cor- 
responds to QCP), 1.0 (Kondo compound) that spans the 
whole range of the solid solution with varied properties. 
In Fig. 2(a), we show x-ray photoemission (XP; hv ~ 
1486.6 eV) valence band spectra collected at 5 K for the 
above four compositions (energy resolution used for these 
spectra = 0.6 eV). The spectra primarily consist of two 
groups of features. The signals corresponding to Ce 6^3/2 
and 5pi/2 photoemission appear at 17.5 eV and 20.3 eV 
binding energies, respectively. The features close to the 
Fermi energy, tp constitute the valence band. Since Ce 



concentration is same in all the samples, the spectra have 
been normalized by the Ce 5p peak intensity. It is ev- 
ident that the spectral weight around 2.8 eV gradually 
reduces with the decrease in Rh content. 

To probe this region with greater clarity, we have ob- 
tained the high resolution (resolution = 0.3 eV) XP spec- 
tra of the valence band region close to as shown in 
Fig. 2(b) after integral background subtraction - a typ- 
ical background function is shown in Fig. 2(a) by solid 
line. The spectral intensity within 1 eV oi tp does not 
change significantly with a; - a small increase is observed 
around 0.5 eV while there is a large decrease in intensity 
around 2.8 eV with the increase in Co content. Such an 
anomalous change can be attributed to the large photoe- 
mission cross-section of Rh 4(i states ('^ 0.012) compared 
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that for Co 3d states (~ 0.0037) In order to visual- 
ize the spectral evolution due to substitution better, we 
renormalized the same spectra by the intensity at 1.5 eV 
binding energy as shown in Fig. 2(c). The increase in 
intensity corresponding to Co 3d photoemission is shown 
by gray shade and the decrease Rh Ad intensity is shown 
by blue shade - the gradual change in spectral intensities 
with the increase in x is beautifully manifested in the fig- 
ure. These results suggest dominant Co 3d contributions 
close to ep, whereas Rh Ad contribute at higher binding 
energies. 

The energy band structure calculations were performed 
for the end members Ce2RhSi3 and Ce2CoSi3 using 
full potential linearized augmented plane wave method 
(FLAPW) within the local density approximations using 
Wien2k softwareJ^ The lattice constants and all other 
structural details were obtained from the literature.— The 
calculated Ce 5p, Rh Ad/Co 3d and Si 3p partial density 
of states (PDOS) for Ce2RhSi3 and Ce2CoSi3 are shown 
in Fig. 2(d) and Fig. 2(e), respectively. The results re- 
veal dominance of the contributions of Rh/Co d states 
in the valence band. The Rh Ad PDOS appear around 3 
eV binding energies consistent with the experimental ob- 
servations. The Co 3d states peak around 1 eV binding 
energy. Clearly, the calculated results correspond well 
with the experimental spectra shown in Fig. 2(c). 

The d PDOS seem to have influenced the conduction 
electron-Ce 4/ hybridization (c/-hybridization) signifi- 
cantly - Ce 4/ PDOS width is broader in Ce2RhSi3 rel- 
ative to that in Ce2CoSi3 as shown in the inset of Fig. 
2(e). Moreover, the intensity of Co 3d PDOS in the Ce 
4/ region is stronger than the Rh 4c? contributions in the 
same energy range. These observations suggest stronger 
Co 3d-Ce Af hybridization than Rh 4d-Ce 4/ hybridiza- 
tion. Evidently, proximity of the d states have strong in- 
fluence on c/-hybridization'^'^ in addition to the influence 
of lattice volume change often referred in these systemsjS 

In Fig. 3, we compare the photoemission spectra col- 
lected using He la, He iia and Al Ka photon energies 
at 5 K. All the spectra are normalized by the intensities 
between 1-2.5 eV binding energy range. The features 
in He ii and XP spectra are quite similar except the fact 
that the XP spectra are broader due to higher resolution 
broadening. Two distinct and sharp features appear in 
the vicinity of in the He ii spectra of both Ce2RhSi3 
and Ce2CoSi3i^i^ The He i spectra, however, do not ex- 
hibit these intensities. Considering contributions of var- 
ious electronic states near ep depicted in Figs. 2(d) and 
2(e), and the dependence of their photoemission cross 
section, the sharp features near the Fermi level can be 
attributed primarily to the photoemission from Ce Af 
electronic states. The feature near ep corresponds to 
4/5 /2 photoemission and is a manifestation of the Kondo 
resonance features. The corresponding spin-orbit split 
component, the 4/7/2 final state appears at a relatively 
higher binding energies 260 meV). In Fig. 3(c), we 
show the He 11 spectra of all the compositions at 5 K. 
Although the Rh 4c? intensity gradually decreases with 
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FIG. 3: Valence band spectra of (a) Ce2RhSi3 and (b) 
Ce2CoSi3 at 5 K obtained using different photon energies ex- 
hibiting results consistent with the band structure results, (c) 
He II spectra of Ce2Rhi_a;Coa;Si3 at 5 K for various values of 
X. The high energy resolution employed in this technique 
reveal distinct Kondo features near the Fermi level. 

the increase in x as also observed in the XP spectra, the 
features near ep seems similar in the whole composition 
range. 

Considering significant enhancement of the Ce Af con- 
tributions in the He 11 spectra relative to the other con- 
tributions in the same energy region, one can extract the 
Ce Af spectral function by subtracting the He i spectra 
from the He 11 spectra as often done in similar Ce-based 
systems. The peak of the Kondo resonance feature usu- 
ally appears just above the Fermi energy (unoccupied) 
and the photoemission spectroscopy probes the higher 
binding energy tail of the Kondo peakJ^i^S Thus, the 
extracted 4/ spectral weight by the difference (He 11 - 
He i) depicting the occupied part of the Ce Af PDOS 
would provide a good estimation of the Kondo resonance 
feature. 

The (He 11 - He i) spectra at different temperatures are 
shown in Fig. 4(a) for different compositions separately. 
The spectra of all the compositions exhibit a sharp fea- 
ture near ep - the intensity of the feature reduces gradu- 
ally with the increase in temperature suggesting presence 
of Kondo compensation effect in every case despite the 
fact that the magnetism in these materials are widely 
varied. The intensity of the Kondo resonance feature at 
a particular temperature gradually increases with the in- 
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FIG. 4: (a) The evolution of the extracted Ce 4/ spectral 
functions with temperature for different compositions, (b) 
The spectral weight transfer obtained by subtracting the room 
temperature spectra from the spectra at different tempera- 
tures. 



crease in x. Since the Kondo temperature increases with 
the increase in x, the above observation suggests that 
the intensity of the Kondo resonance feature is directly 
related to the Kondo temperature, Tk-^ While such a 
scenario is expected in a Kondo system, the similar ob- 
servation in long range ordered systems (x < 0.6) indi- 
cates the existence of a finite Kondo temperature scale 
even in these cases. 

The Kondo coherence temperature for the intermedi- 
ate compositions estimated^ from the electrical transport 
properties is ^ 8 K for x = 0.3 and ~ 50 K for x — 0.6. 
When the measurement temperature is below the coher- 
ence temperature, there will be additional complications 
due to the expansion of the Fermi volume and presum- 
ably manifested in the evolution of the Kondo resonance 
feature. The Kondo temperatures for the compositions, x 
= 0.0, 0.3 and 0.6 are 12 K, 30 K and 40 K, respectively,^ 
The photoemission spectra collected at 120 K and 150 K 
are well above the coherence temperature, Kondo tem- 
perature and Neel temperature for all the compositions 
and still exhibit finite intensity of the Kondo resonance 
feature. At all the temperatures, we find qualitatively 
similar trend of the evolution of the Kondo resonance 
feature across the solid solution indicating that the spec- 
tral weight redistributions, if any, due to the onset of the 
Kondo coherence are not altering the inferences drawn. 
This manifests applicability of the single impurity Kondo 



physics in the electronic structure^ and the Kondo tem- 
perature scale remains finite at the QCP. 

The intensity of the corresponding spin-orbit split fea- 
ture (4/7/2 signal) appearing around 260 meV binding 
energy (marked by a vertical dashed line in Fig. 4(a)) 
in each of the samples also reduces as the temperature is 
increased. Interestingly, the intensity of the 4/7/2 signal 
seems to be similar in all the compositions although the 
intensity of the Kondo resonance feature at the Fermi 
level increases with the increase in Co concentration. At 
finite temperature, the spectral intensity representing the 
occupied part of the density of states becomes lower at 
due to the thermal effect represented by the Fermi-Dirac 
distribution function. The Fermi-Dirac distribution func- 
tion depletes the intensity below ep and enhances the in- 
tensity above Ci;- in a symmetric manner with the increase 
in temperature. We have subtracted the room temper- 
ature spectra from the spectra at different temperatures 
(/(e,T) — /(e,300ii')) to investigate the spectral weight 
transfer with temperature. Such exercise would provide 
a change in spectral weight symmetric with respect to 
the Fermi level in a Fermi liquid system. The subtracted 
spectral functions, shown in Fig. 4(b), however, exhibit 
asymmetric spectral evolutions with temperature - the 
gain of spectral weight below the Fermi level is larger 
than the loss of weight above the Fermi level. Such an 
effect appear in the presence of a peak such as Kondo 
resonance peak in the vicinity of the Fermi level. The 
asymmetry seems to be weaker towards Rh end indicat- 
ing weaker Kondo peak. This clearly indicates that the 
temperature dependence of the Kondo peak and the spin- 
orbit peak is different. 

To investigate the spectral changes with composition, 
a; at a particular temperature, we re-plot the Ce 4/ spec- 
tral functions as a function of composition at a particular 
temperature in Fig. 5(a). The experimental spectra ex- 
hibit gradual increase in intensity of both the peaks with 
the increase in x that can be associated to the change 
in the Kondo temperature as discussed above and also 
observed earlier in other Ce-based compounds,^ 

We compare the evolution of the intensities of the 
Kondo resonance peak and its spin-orbit split compo- 
nent in Fig. 6. In order to deconvolve the Fermi Dirac 
distribution induced depletion of the spectral density of 
states at the Fermi level, we symmetrized the spectra 
with respect to the Fermi level - if /(e) represents the 
spectral function, the symmetrized spectra can be ex- 
pressed as /(e) = /(e — ep) + I{eF — e)- It is often found 
that the intensity at the Fermi level can be estimated 
quite accurately by such symmetrization process, which 
is independent of the presence/absence of particle-hole 
symmetry. The symmetrized spectra are shown in Fig. 
5(b) those represent the intensity of the Kondo resonance 
feature (KRF). 

The estimated intensities of the spin orbit peak (SOC) 
obtained from Fig. 4(a) and the Kondo peak (KRF) 
obtained from Fig. 5(b) are shown in Fig. 6(a) 
(for CeaRhSia), 6(b) (for Ce2Rho.7Coo.3Si3, 6(c) (for 
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FIG. 5: (a) The evolution of the extracted Ce 4/ spectral 
functions with composition, x at different temperatures, (b) 
The spectral functions obtained by symmetrizing the experi- 
mental spectra to estimate the spectral weight the Fermi level. 

Ce2Rho.4Coo.6Si3 and 6(d) (for Ce2CoSi3). Interestingly, 
the compositions exhibiting long range order show weaker 
temperature dependence of the Kondo resonance feature 
(KRF) compared to the spin orbit coupled (SOC) peak 
- the trend gets reversed on the other side of the phase 
diagram. This is shown by plotting the relative inten- 
sity ratio (KRF/SOC) in the insets of Fig. 6(e) in the 
corresponding phase regime. In the weak hybridization 
limit, the intensity ratio, KRF/SOC decreases with the 
decrease in temperature - the Kondo peak grows less 
rapidly than the spin orbit peak. On the other side 
(strong coupling limit), the ratio, KRF/SOC gradually 
increases with the decrease in temperature. Interestingly, 
the trend at the quantum critical point follows the strong 
coupling behavior. 

Kondo effect essentially arises from the hybridization 
of the conduction electrons with the impurity moment. It 
is well known that in the strong coupling limit, the elec- 
trons corresponding to the local moment forms a quan- 
tum mechanically entangled states with the conduction 
electrons, termed as Kondo singlet. In such a scenario, 
the Fermi surface expands due to the additional contri- 
bution from the electrons forming the local moment - 
here the 4/ electrons. The Kondo temperature essen- 
tially provides an energy scale deriving the expansion of 
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FIG. 6: Intensity variation of the Kondo resonance feature 
(KRF) and its spin orbit coupled (SOC) component as a func- 
tion of temperature of (a) Ce2RhSi3, (b) Ce2Rho.7Coo.3Si3, 
(c) Ce2Rho.4Coo.6Si3 and (d) Ce2CoSi3. (e) A schematic 
phase diagram showing the long range ordered region (yel- 
low) and Kondo compensated region separated by a quantum 
critical point (QCP). The left and right insets show the in- 
tensity ratio of KRF/SOC as a function of temperature for 
corresponding compositions in Ce2Rhi_xCoa;Si3. Tjv and T* 
represent the Neel temperature and Kondo coherence temper- 
ature and are obtained from Ref. [6]. 



the Fermi volume.— In many cases, it is observed that 
the Kondo temperature scale vanishes at the approach to 
QCF.^"^'^^ The results in the present study reveal finite 
Kondo temperature scale at the QCP in the solid solu- 
tion Ce2Rhi_j;Coa;Si3 suggesting that the magnetic Ce 
4/ electrons near QCP acquire itineracy at lower tem- 
peratures - the antiferromagnetism near the QCP occurs 
among the itinerant electrons and hence indicates the ap- 
plicability of an SDW picture in these cases.=S'i& Interest- 
ingly, the results in Fig. 6(e) suggests that if the Kondo 
temperature is higher than Neel temperature, the in- 
tensity of the Kondo resonance feature exhibits stronger 
temperature dependence compared to its spin-orbit split 
counterpart. The scenario gets reversed, if the Kondo 
temperature is lower than the ordering temperature. 



6 



IV. CONCLUSIONS 

In summary, we studied the evolution of the electronic 
structure of Cc2Rhi_.j,Coj.Si3 using high resolution pho- 
toemission spectroscopy. The signature of Kondo reso- 
nance feature is revealed in the high resolution low tem- 
perature spectra of all the compounds providing evidence 
of finite Kondo temperature scale at the quantum criti- 
cal point. Interestingly, the temperatiu-c dependence of 
Kondo resonance feature is found to be stronger than 
its spin-orbit split counterpart if the Kondo temperature 
is higher than the magnetic ordering temperature. The 
spin-orbit part show stronger temperature dependence if 
the Kondo temperature becomes lower than the ordering 
temperature. While finite Kondo temperature scale in 



this regime establishes applicability of the spin density 
wave picture at the approach to the quantum critical- 
ity, the spin-orbit coupling appears to play a significant 
role in this regime. These results provide an insight to 
the origin of quantum fiuctuations often observed in the 
proximity of the onset of exotic phases like superconduc- 
tivity. 



V. ACKNOWLEDGEMENTS 

One of the authors S.P., thanks the Council of Scien- 
tific and Industrial Research, Government of India for 
financial support. 



* Corresponding author: kbmaiti@tifr.rcs.in 

1 J. A. Hertz, Phys. Rev. B 14, 7183 (1976). 

2 S. Doniach, Physica B & C 91, 231 (1977). 

^ Q. Si, S. Rabello, K. Ingersent, and J. L. Smith, Nature 

413, 804 (2001). 
" A. J. Millis, Phys. Rev. B 48, 7183 (1993). 
^ M. Garnier, K. Breuer, D. Purdie, M. Hengsberger, Y. 

Baer, and B. Dcllcy, Phys. Rev. Lett. 78, 4127 (1997). 
® S. Patil, Kartik K. Iyer, K. Maiti, and E. V. Sampatliku- 

maran, Phys. Rev. B 77, 094443 (2008). 

R. A. Gordon, C. J. Warren, M. G. Alexander, F. J. DiS- 

alvo, and R. Pottgen, J. Alloys Compds. 248, 24 (1997). 
® S. Majumdar, M. M. Kumar, R. Mallik, and E. V. Sam- 

pathkumaran. Solid State Comm. 110, 509 (1999). 
^ S. Patil, V. R. R. Medicherla, R. S. Singh, S. K. Pandey 

E. V. Sampathkumaran, and K. Maiti, Phys. Rev. B 82, 

104428 (2010). 

1° S. Patil, S. K. Pandey, V. R. R. Medicherla, R. S. Singh, 
R. Bindu, E. V. Sampathkumaran and K. Maiti, J. Phys.: 
Condens. Matter 22, 255602 (2010). 

" K. Maiti, U. Manju, S. Ray, P. Mahadevan, l.H. Inouc, 
C. Carbonc, and D.D. Sarma, Phys. Rev. B 73, 052508 

(2006) ; K. Maiti, A. Kumar, D.D. Sarma, E. Weschke, 
and G. Kaindl, Phys. Rev. B 70, 195112 (2004); K. Maiti 
and D. D. Sarma, Phys. Rev. B 61, 2525 (2000). 

^2 K. Maiti and R.S. Singh, Phys. Rev. B 71 161102(R) 
(2005); K. Maiti, R.S. Singh, and V.R.R. Medicherla, 
Phys. Rev. B 76, 165128 (2007); R. S. Singh, V. R. R. 
Medicherla, and K. Maiti, Appl. Phys. Lett. 91, 132503 

(2007) . 

" T. Nakano, K. Sengupta, S. Rayaprol, M. Hedo, Y. Uwsr 
toko, and E. V. Sampathkumaran, J. Phys.: Condens. 
Matter 19, 326205 (2007). 

I. Das and E. V. Sampathkumaran, J. Magn. Magn. Mater. 
137, L239 (1994). 

K. Mukherjee, K. K. Iyer, S. Patil, K. Maiti, and E.V. Sam- 
pathkumaran, J. Phys. Conf. Series 273, 012010 (2011); K. 



Maiti, S. Patil, G. Adhikary, and G. Balakrishnan, J. Phys. 
Conf. Series 273, 012402 (2011); K. Maiti, S. Patil, and E. 
V. Sampathkumaran, AIP Conf. Proc. 1347, 181 (2011) 
J. M. Lawrence, P. S. Riseborough, and R. D. Parks, Rep. 
Prog. Phys., 44 1 (1981). 

P. Blaha, K. Schwarz, G.K.H. Madsen, D. Kvasnicka, and 
J. Luitz, WIEN2k, An Augmented Plane Wave -I- Local 
Orbitals Program for Calculating Crystal Properties (Karl- 
heinz Schwarz, Tcchn. Universitat Wien, Austria), 2001. 
ISBN 3-9501031-1-2. 

J.J. Yeh and I. Lindau, At. Data Nucl. Data Tables 32, 1 
(1985). 

F. Patthey, W. -D. Schneider, Y. Bacr, and B. Dcllcy 
Phys. Rev. Lett. 58, 2810 (1987); F. Patthey, J. M. Imer, 
W. -D. Schneider, H. Beck, Y. Baer, and B. Delley, Phys. 
Rev. B 42, 8864 (1990). 
^° D. Ehm, S. Hiifner, F. Reinert, J. Kroha, P. Wolfle, O. 
Stockert, C. Geibel, and H. v.Lohneysen, Phys. Rev. B 76, 
045117 (2007); M. Klein, A. Nuber, F. Reinert, J. Kroha, 
O. Stockert, and H. v.Lohneysen, Phys. Rev. Lett. 101, 
266404 (2008). 

O. Gunnarsson and K. Schonhammer, Phys. Rev. B 28, 
4315 (1983). 

O. Gunnarsson, N. E. Christensen, and O. K. Andersen, J. 
Magn. Magn. Mater. 76-77, 30 (1988); doi: 10. 1016/0304- 

8853(88)90309-5. 

M. Klein, A. Nuber, F. Reinert, J. Kroha, O. Stockert, and 
H. V. Lohneysen, Phys. Rev. Lett. 101, 266404 (2008). 
S.-I. Fujimori, Y. Saitoh, T. Okane, A. Fujimori, H. Yam- 
agami, Y. Haga. E. Yamamoto, and Y. Onuki, Nat. Phys. 

(London) 3, 618 (2007). 

P. Gegenwart, Q. Si, and F. StegUch, Nat. Phys. (London), 
4, 186 (2008). 

^® P. Gegenwart, T. Westcrkamp, C. Krcllner, Y. Tokiwa, S. 
Paschen, C. Geibel, F. Steglich, E. Abrahams, and Q. Si, 
Science, 315, 969 (2007). 



